function EQ = eqsystem(A, B, C, s)
EQ1= ((A.'*B)^(-1))*diag(A*B.')-C;
ini = ones(s,1);
EQ2=(ini.'*B)-1;
EQ = [EQ1; EQ2];



